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CN i Abstract 



In this article, we analyze a two-level finite element method for the equations of motion arising 
in the flow of 2D Oldroyd model with non-smooth initial data. It involves solving the non-linear 
problem on a coarse grid of mesh-size H and solving a linearized problem on a fine grid of mesh- 
CNJ 1 size h, h << H. The method gives optimal convergence rate for velocity in ii^-norm and for 

pressure in L 2 -norm. The analysis takes in to account the loss of regularity of the solution of the 
Oldroyd model at initial time. 
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We consider a two-level semi-discrete Galerkin approximations to the following system of equations 
l/-^ , of motion arising in the Oldroyd fluids of order one (see [26J ) : 



(1.1) ^ + u-Vu-/iAu- / P(t-T)Au(x,T)dT + Vp = f(x,t), x e n, t > 
&t Jo 

with incompressibility condition 

(1.2) V • u = 0, x G O, t > 0, 



X 1 

5^ " and initial and boundary conditions 



(1.3) u(x,0) = uo in f2, u = 0, on dQ, t>0. 

Here, is a bounded domain in M 2 with boundary dQ, fi = 2k\~ 1 > and the kernel j3{t) = 
7exp(— 6t), where 7 = 2\~ 1 (v — kX^ 1 ) > and S = A -1 > 0. u and p are the velocity field and 
pressure, respectively. Further, the forcing term f and the initial velocity uo are given functions in 
their respective domains of definition. For more details, we refer to [8] and references cited therein. 

There is a considerable amount of literature devoted to Oldroyd model by Oskolkov, Kotsiolis, 
Karzeeva, Sobolevskii etc, see [31 HI [TH [HI [28] and recently by Lin, He et al. [HI [HI [33] and by 
Pani et al. [291 SO] ■ A brief introduction on the continuous and semi-discrete cases can be found in 
[8]. In fact, in [8], we have established a priori estimates and regularity results for the solution pair 
{u,j?} of (jl.ip - (jl.3p under realistically assumed initial data and when f, f t G L°°(1R + ; L 2 (il)). 

Semi-discrete Galerkin finite element approximations and error analysis for this problem are 
carried out in [HI [291 E] , while fully discrete cases can be found in [91 [301 E3] • Other works on this 
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problem, that is, (jl.ip - (jl.3p . can be attributed to He, Wang et al., see [31] [32] [33] . But there is no 
work on two-level or similar methods for this problem, to the best of author's knowledge and this is 
the first work in this direction. 

Two-level or two-grid methods are well-established and efficient methods for solving non-linear 
partial differential equations. It involves solving the non-linear problem on a coarse grid (rather 
than on a fine grid, thereby computationally cost effective), and solving a linearized problem on a 
fine grid. In other words, in the first step, we discretize the non-linear PDE on a coarse mesh, of 
mesh-size H and compute an approximate solution, say, u#. Then, in the second step, we formulate 
a linearized problem, out of the original one, using u# and discretize it on a fine mesh, of mesh-size 

h, h << H, thereby, compute an approximate solution, say, u h . With appropriate h,H, we obtain 
same order of convergence of the error u — u h , as that of the error obtained by semi-discrete finite 
element Galerkin approximation on the find grid; but with far less computational cost, since, instead 
of solving a large non-linear system, we solve a small non-linear system and a large linear system. 

In this article, we study the following two-level finite element approximation for the problem 
(jl.ip - (|1.3p : First, we compute a semi-discrete Galerkin finite element approximations (uh,Ph), over 
a coarse mesh of mesh-size H. Then, we use the approximation u# to compute a semi-discrete 
Galerkin finite element approximations (u h ,p h ) of the following linearized Oldroyd problem: 

(1.4) Ui-i/Au + u/fVu/f- / j3(t - s)Au(s) ds + Vp = f , V • u = 0, in U 

Jo 

over a fine mesh of mesh-size h, h « H. 

The above algorithm is nothing new and in fact, this and similar algorithms have been studied on 
numerous occasions for both Stokes and Navier-Stokes problems. But to the best of our knowledge, 
no study has been done for our problem and with non-smooth initial data (even for Navier-Stokes), 

i. e., initial velocity is in Hj and not in Hq n H 2 or in higher order Sobolev space. Our main aim in 
this work is to do error analysis of this two- level method, under non-smooth initial data. 

Two-grid method was first introduced by Xu |35[ [36] for semi-linear elliptic problems and by 
Layton et al. [20], [2T] [22] for steady Navier-Stokes equations. It was carried out for time dependent 
Navier-Stokes by Girault and Lions [7] for semi-discrete case. The method may vary depending on 
the linearize problem, to be solved in the second step; like, in the case of Navier-Stokes, one can 
chose a Stokes problem or an Oseen problem or a Newton step to solve on the fine mesh. Several 
works in this direction, involving both semi-discrete and fully discrete analysis, can be found in 
[D[2l[IIl[l2l[l5l[23l[Ml[25l[27]and references therein. 

We would like to note here that, similar to Navier-Stokes (see |16] for a discussion), the regularity 
of the solution of our problem breaks down as t — > 0. To assume otherwise, the initial data have 
to satisfy certain non-local compatibility conditions that are very difficult to verify and physically 
irrelevant. We have avoided these conditions in jS], while carrying out semi-discrete error analysis, 
thereby implying breakdown of regularity of the solution at initial time. Here also, we have worked 
out under similar circumstances. 

Recently, two articles [5] [6], in the context of Navier-Stokes, have taken this "loss of regularity 
of the solution at initial time" into account. For example, both these articles assume no more than 
second-order spatial derivative of the velocity bounded in L 2 , up to initial time t = 0. But in our 
case, due to non-smooth initial data, we have ||u(t)||2 ~ 0(i _1//2 ), see [8] Theorem 3.1]. In other 
words, no more than first-order spatial derivative of the velocity is bounded in L 2 , up to initial time. 
In this article, we present the following results: 

\\u h (t) - u h (t)\\i + \\ Ph (t)-p h (t)\\ < Kt-'H 2 , 

where (u h ,p h ) are the solutions of the two- level method on a grid of mesh-size h and (u/j,^) are 
the solution of the standard Galerkin finite element, which satisfy the following: (see [8]) 

||u(t) - u h (t)||i + \\p(t)- Ph (t)\\ < Kt-^h. 
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Therefore, away from 0, this two-level method produces result similar to standard Galerkin finite 
element for H = 0{h 1 / 2 ). This implies that we obtain similar result at a less computational cost, 
since this two-level method treats the non-linearity only in a coarse grid rather than a fine grid 
treatment as in Galerkin finite element. 

Apart from the error analysis, we have also looked into the stability of the method, by working 
out the error analysis. We would like to note here that a proper justification can only be done with 
the analysis of full discretization and numerical implementation. This will be done in our next effort. 

In the following sections, we will assume C and K to be generic positive constants, where K 
generally depends on the given data, that is, uo,f. And for the sake of convenience, we would write 
v(t) as v when there arises no confusion. 

The article is organized as follows. In section 2, we introduce some notations and preliminaries. 
Section 3 deals with Galerkin finite element approximations, whereas two- level method is discussed 
in Section 4. Error analysis for two-level is carried out in Section 5. 

2 Preliminaries 

For our subsequent use, we denote by bold face letters the M 2 -valued function space such as 

Hi = [iT W, L 2 = [L 2 (fl)] 2 . 
Note that Hj is equipped with a norm 

2 2 

||Vv|| = ( £(^,^)) 1/2 = (^2(Vv u Vv t )) 1/2 . 

i,j=l i=l 

Further, we introduce the following divergence free function space: 

Ji = {</> € Hq : V • = 0}. 

For any Banach space X, let L p (0, T; X) denote the space of measurable X -valued functions <f> on 
(0, T) such that 

/ \\cj>{t)\\ p x dt < oo if 1 <p < oo, 
Jo 

and for p = oo 

ess sup < oo if p = oo. 

0<t<T 

Through out this paper, we make the following assumptions: 

(Al). For g G L 2 , let the unique pair of solutions {v <G 3\,q € L 2 /R} for the steady state Stokes 
problem 

-Av + Vq = g, V • v = in Q; v\ dn = 0, 
satisfy the following regularity result 

IMk + WqWhi/r < C||g||- 

(A2). The initial velocity uo and the external force f satisfy for some positive constant Mo, and for 
T with < T < 00 

u <E Ji, f,f t ,f« G L°°(0,T;L 2 ) with ||u ||i < M , sup {||f ||, ||f t ||, ||f«||} < M . 

0<t<T 
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Before going into the details, let us introduce the weak formulation of (jl.ip - (jl.3p . Find a pair of 
functions {u(f),p(i)}, t > 0, such that 

(u t , 0) + /i(Vu, V0) + (u • Vu, 0) + f p(t - s)(Vu(s), V0) = (p, V ■ 0) + (f , 0) V0 G Hj, 

jo 

(2.1) (V-u, X )=0 V X GL 2 . 
Equivalently, find u(i) G Ji, t > such that 

(2.2) (u t , 0) + ^(Vu, V0) + (u • Vu, 0) + / - S )(Vu(s), V0)ds = (f, 0), V0 G Ji. 

For existence and uniqueness, and the regularity of the solution of the problem (|2.ip or (|2.2p . 
we refer to [8]. In fact, in our subsequent section, we will use various estimates of u and p, without 
explicit mention of them. These estimates can be found in Lemma 3.1-Theorem 3.2, [SJ. 
We present below the positivity of the kernel /3, which will be used in our subsequent analysis. The 
result is borrowed from [8j Lemma 2.1]. 

Lemma 2.1. For arbitrary a > 0, t* > and 4> G L 2 (0, t*), the following positive definite property 
holds 

exp [-a(t - s)]<f>(s) ds ] <f>(t) dt > 0. 



Lemma 2.2 (Gronwall's Lemma). Let g,h,y be three locally integrable non-negative functions on 
the time interval [0, oo) such that for all t > 

y(t) + G(t) < C + f h(s) ds+ [ g(s)y(s) ds, 
Jo Jo 

where G(t) is a non-negative function on [0, oo) and C > is a constant. Then, 

y(t) + G{t) <(C + J h{s) dsjexp^J g(s) da). 

3 Galerkin Finite Element Method 

From now on, we denote h, with < h < 1, to be a real positive spatial discretization parameter, 
tending to zero. Let and be two family of finite element spaces, finite dimensional subspaces 
of Hj and L 2 /R, respectively, approximating velocity vector and the pressure. Assume that the 
following approximation properties are satisfied for the spaces and L^: 

(Bl) For each w G Hq n H 2 and q G H 1 /M, there exist approximations i^w G H/j and ]\q G such 
that 

||w - i h w\\ + h\\S/(w - i h w)\\ < C7i 2 ||w|| 2 , \\q - j h q\\ < Ch\\q\\i. 
Further, suppose that the following inverse hypothesis holds for w/j G H^: 

(3.1) ||Vw A || < Ch^Wwhl 

For v, w, <p G Hq, set 

o(v, 0) = (Vv, V0), 6(v, w, 0) = -(v • Vw, 0) - -(v • V0, w). 

Note that the operator &(•, •, •) preserves the antisymmetric property of the original nonlinear term, 
that is, 

b(v h , w/j, w h ) = Vv^, w h G H h . 
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The discrete analogue of the weak formulation (|2.ip now reads as: Find Uh(t) G and ph(t) G Lh 
such that u?j(0) = uoh and for t > 

(ufct> </>/») + f^a(u h ,4> h ) + b(u h ,u h ,4> h )+ [ fi(t - s)a(u h (s),4> h ) ds - (p h ,V ■ 4> h ) = (f,(f> h ), 

Jo 

(3.2) (V.u h ,xh) = 0, 

for <p h G H/j, x/j € Lft. Here uo/j G H/j is a suitable approximation of uo G Ji. Or we can define 

3h = {"h G : (xh, V ■ = Vx/, G 4}, 
for an equivalent weak formulation: Find u/ l (t) G J/i such that u/ l (0) = uoh and for t > 

(3.3) (u M ,4>h) + (Ja(uh,<f> h )+ P(t-s)a(uh(s),<f> h ) ds = -b(u h ,u h ,<ff h ) + (f,<f> h ) V0 h G J ft . 

J o 

For the well-posedness of (|3.2[) or (|3.3p . we refer to [HJ [29]. Note that the pressure Ph{t) G is 
unique up to a constant and unique in a quotient space L^/N^, where 

N h = {q h G L h : {q h , V • fc ) = for 0^ G HJ. 

The norm of L^/N^ is given by 

hhhyN h = inf ||gh + X/i||- 

For continuous dependence of the discrete pressure Ph(t) G L^/Nh on the discrete velocity Uh(t) G J^, 
we assume the following discrete inf-sup (LBB) condition: 

(B2') For every G Lh, there exists a non-trivial function cf) h G such that 

Kfffci V • > CIlV^fcllHghllia/jVfc, 
where the constant C > is independent of fa. 

Moreover, we also assume that the following approximation property holds true for J/,. 
(B2) For every w G Ji H H 2 , there exists an approximation r/,w G Jh such that 

||w — r/jw|| + /i||V(w — r/jw)|| < C/i 2 ||w||2. 

The L 2 projection Ph : L 2 i-> J/j satisfies the following properties (see [IS]): for 6 J/,, 

(3.4) ||0 - f\0|| + /i||VP h 0|| < C/i||V0||, 
and for G Ji n H 2 , 

(3.5) ||0 - P h (f>\\ + V(0 - P h <f>)\\ < Cfa 2 ||A0||. 

We now define the discrete operator A^ : i— » through the bilinear form a(-, •) as 

(3.6) a(v h ,4>h) = (-A h v h ,<j>) Vv h ,4> h G H ft . 

Define = P/ l (— A^). The restriction of A/j to J/i is invertible and we denote the inverse by A^ 1 . 
Since — A^ is self-adjoint and positive definite, we define discrete Sobolev norms on Jh as follows: 



v h\ 



-A ft ) r / 2 v ft ||, v h G Jh, T G 



We note that, in particular, ||v^||o = \\vh\\ an d ll v ftlli = l|V v /ill f° r v h G Jh, and || • H2 and ||A^ • 

are equivalent on Jh- For further detail, we refer to [I6| [T7], 

Below, we present some estimates of the non-linear term (see (3.7) from |17J). 



5 



Lemma 3.1. Suppose conditions (A.1), (Bl) and (B2) are satisfied. Then there exists a positive 
constant K such that for v, w, 4> £ H^, the following holds: 



(3.7) 



|(v Vw,0)| <C< 



V2|| Vv ||l/2|| Vw ||l/2|| AhW ||l/2|| || ) 
V2 1/3 1| VW|| , 



1 / 2 ||Vv|] 1 /2||Vw|||| 1 
llVwIIII^II^HA^ 
llVwf/ailAfcwII 1 ^ 



6|| 1 /2|| V 0|| 1 /2 ) 
11/2 



|l/2 



l|V0|| 



1/2 



Examples of subspaces H/j and satisfying assumptions (Bl) and (B2') are abundant in literature 
(see [8] and the references cited there). We present below a lemma dealing with a priori and regularity 
estimates of u^. 

Lemma 3.2. Let the assumptions of Lemma \3.1\ hold. Additionally, let < a < min{^Ai/2, 5} The 
semi- discrete Galerkin approximation u/j of the velocity u satisfies, for t > 0, 



(3.8) 



2 + e -2at 



e 2as \\u h (s)\\i ds < K, 



(3-9) \\u h (t)\\l + e- 2at / e 2as {\\u h (s)\\ 2 + \\u Ks (s)\\ 2 } ds < K, 

Jo 

(3.10) (r*(t)) 1/2 {K(t)|| 2 + \\ Ph (t)\\ H i /R } < K, 

where T*(t) = min{l,t} and K depends only on the given data. In fact, K is independent of h. 

Remark 3.1. Here \\ > is the least eigenvalue of the Stokes operator (see JBj). 

Proof. The estimates can be found in Lemma 4.2 of [8]. The additional estimates of u^t and of p 
can be derived as in the proof of Theorem 3.1 of [8]. □ 

Lemma 3.3. Under the assumptions of Lemma \3.S\ the semi-discrete Galerkin approximation 
of the velocity u satisfies, for t > and for r € {0, 1}, i E {1,2}, r + i < 2, 



(3.H) 



r+2i-li 



D\u h {t)f r + e 



-2at 



r+2i-12as 



where D\ = And 



-2at 



2„2as i 



(s)\\ 2 + 



e Zas \\D l s u h (s)\\ z r+1 ds<K, 



(3.12) e ~" / {T~[s)re"-{\\u ht . 

Jo 

(3-13) (r*(t)) 3/2 {||u M || 2 + \\p h ,t\\ m/R } < K. 

Here K depends only on the given data. In particular, K is independent of h and t. 

Proof. The proofs of (|3.1ip for i = 1 and of (13.12j) follow in a similar fashion as that of Theorem 
3.2 from [8]. For the remaining estimates, we need to differentiate (|3.3|) twice with respect to time 
and choose cj) h = (r*(t)) 3 e 2a *u/ litt and again follow the ideas of the proof of Theorem 3.2 from [5J to 
complete the rest of the proof. □ 

The following semi-discrete error estimates are proved in 



Theorems 5.1 and 6.1] . 



Theorem 3.1. Let Q, be a convex polygon and let the conditions (A.\)-(A2) and (fS\)-(SS2) be 

satisfied. Further, let the discrete initial velocity uoh € Jh with uq^ = -f\uo, where uq 6 Ji. Then, 
there exists a positive constant K , that depends only on the given data and the domain Vt, such that 
for < T < oo with t G (0, T) 

||(u - u h )(t)|| + h\\V{M - u h ){t)\\ + h\\(p- Ph )(t)\\ < Ke Kt h 2 t- l l 2 . 
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For the error analysis of the two-level method, we need a couple of additional estimates of the 
error due to Galerkin finite element, that is, of (u — Uh)(t). In fact, we need estimates of the 
time derivative of the error. As in the case of parabolic problem, we need to differentiate the error 
equation, that is, the equation in (u — Uh)(t), with respect to time and derive several estimates to 
conclude the result. It is technical and lengthy, but follow the proof techniques of Section 5 from [8j. 
For the sake of completeness, we try to sketch a proof here. We will be as brief as possible, conveying 
only the ideas, just to keep it clean and concise. 

Lemma 3.4. Under the assumptions of Theorem \3.1[ the following result holds: 

(3.14) e~ 2at f<7i( a )[|(u.-u hja )( a )|| 2 (fc < K(t)h\ 

Jo 

where <7i(t) = (r* (i)) 2 e 2at and K(t) is in fact, of the form Ke , but is independent of h. a is given 
by Lemma \3.2l 

Proof. We first recall the equation of E = u — (see [El (6.2)]). 

(3.15) (E t ,0 fc ) + /io(E,0 fc )+ / P(t-s)a(E(s),<t> h ) ds = A lih (0 h ) + (p, V • V0 h € 3 h , 

Jo 

where 

(3.16) Ai >h (0 h ) = b(u h , u h , 4> h ) - 6(u, u, <j> h ) = -b(u h , E, <j> h ) - 6(E, u, <j> h ). 
We differentiate (|3.15p to find 

(3.17) (Ett, 4>h) + Va(Vt, <t> h ) + 7 a(E, <f> h ) + / /3 t {t - s)a(E(s), <j> h ) ds = Ai >h , t (0 h ) + (p t , V • h ), 

Jo 

where 

(3.18) A w (0 h ) = -6(u M , E, h ) - 6(u Al E t , h ) - 6(E t , u, - 6(E, u t , <j> h ). 

With as the solution of a linearized Galerkin approximation (5.2) of [8], we break the error Ef in 
linear and non-linear parts, as 

E t = (u - \ h ) t + (y h - u h ) t := £ t + r) t . 

The equation in £ t is given by 

(3.19) (|tt, <t> h ) + $ h ) + 7«(€, + / Pt(t-s)a(€(s),(t>h) <k = (pt, V ■ (t> h ). 

Jo 

With P fe : L 2 (fi) — >■ J/, as L 2 -projection, we choose <f> h = cri(t)P h £ t = ai(t)(£ t — (u — Pftii) t ) in (|3. 19[) 
to find (here a^t) = (t* (t)) 2 e 2at , r*(t) = min{l,t}) 

~{^i(*)ll€*ll 3 } + A"i(*)ll€*ll! = ^l,*(*)ll^f " 7*t(*)««,6 " ( u " fli«)t) 

(3.20) +*l (u " ^u) t ) + A«ri(*)o(e t> (u - Pftu)t) 

-ai(t) / /3 t {t-s)a(Z(s),Z t -(u-P h u) t ) ds + <ri(t)(p t ,V • Pft£ t ). 

JO 
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We note that 



d II/.. n ..\ 112 



<Ti(t)(£tt> ( u " p h")t) = «r x (t)((u - Pfcii)tt, (u - P ft u) t ) = "^^IKu - P^u) 

= |{^ll(u-P,u) t || 2 }-^||(u-P /l u) 4 || 2 , 
ai(t)( Pt , V • P ft £ t ) = (Ti(t)(pt - ifcpt, V • P^ t ). 

Incorporate these in (I3,20p , Then, use projection properties (Bl) and Cauchy-Schwarz inequality to 
have 

~Mt)Ut\\ 3 } + W(t)Ut\\l < Ca(t)Ut\\ 2 + Ca^UhiUth + cfc||u*|| a } 
|{^||(u - P h u) t \\ 2 } - ^||(u - P h u) t \\ 2 + Ch ai(*)||e,||i{||ii«||2 + Mil} 

(3.21) +Cai (t){U t \\ 1 + ch\\u t \\ 2 } [ p t ( t -s)U(s)\\id S . 

Jo 

(a(t) = r*(t)e 2at , r*(t) = min{l,t}) 

Use Young's inequality and kickback argument. Integrate the resulting inequality. We handle the 
double integral term as (4.6) of [29] to find 



<ri(t)Mt\\ 2 + H <n(s)\\Us)\\lds<C a(sMM\\ 2 ds + Ch 2 <n(s){||u fl ( 5 )||i + \\ Ps (s)\\ 2 } ds 
Jo Jo Jo 

c[ e^Uis^jds + Cha^WutWl + Ch [ ||u s ( S ) ||f ds. 
Jo Jo 

Use Theorem 3.2 from [8] to bound the terms involving u and p. We borrow the estimate (5.7) from 
|29j (although the estimate is for f = 0, can be achieved similarly in our case), that is, 

(3.22) e~ 2at f e 2as U{s)\\\ds < Kh 2 

Jo 

to observe that 

(3.23) (r*(t)) 2 U t f + f,e- 2at f a^UM^ds < Ce~ 2at f a(s)UM\\ 2 ds + Kh 2 . 

Jo Jo 

To estimate the first term on the right hand-side of (|3.23p . we first recall the equation in £ (see [HJ 
(5.3)]): 

(£t>0h) +A*a(£,0 ft ) + f s)a(£(s),<f> h ) ds = (p,V • <f> h ), V0 h G 3 h . 

Jo 

Choose 4> h = a{i)Ph£t above. 

<r{t)Utt + ~ !W*)ll£Hi} < Ce 2at U\\ 2 + Ch.a{t){\\i t \\\\*th + Uhh*h + hhUM 

+Ca(t){U t \\ 1 + Ch\\vi t \\ 2 } [* Pt(t-8)M(s)\\i ds. 

Jo 

Integrate with respect to time, use kickback argument, Lemma 13.31 and (|3.22p to conclude 

(3.24) e~ 2at [\(s)MM\\ 2 ds + T*(t)U\\i<Kh 2 . 

Jo 



We have again used Theorem 3.2 from [8j to bound u t and p. Using (|3.24p in (|3.23p leads us to 

(3.25) (r*(t)) 2 U t \\ 2 + e-^ f ' o x (s)\\U°)td8 < Kh 2 . 

Jo 

Next, we use parabolic duality argument (similar to the arguments (5.8)-(5.15) and the rest, in |29j ) 
to establish 

(3.26) e~ 2at f <ri(8)\\Us)\\ 2 d8<Kh A . 

Jo 

An estimate of r\ t would now complete the proof. By definition, we can easily deduce the equation 
satisfied by rj t . 

(3-27) (r} tt , 4> h ) + fia{r] t , tf> h ) + ja(rj, cj> h ) + f f3 t (t - s)a(rj{s), (j> h ) ds = A w (0,J, 

Jo 

where 4>\ h,t ls given by (13.18p . Choose <fi h = ai(t)(—Ah)~ 1 rj t and similar to (I3.2ip . we find 

^^ l(t)ll ^ l| - l} + ' UCTl(t)l|r?il12 - + c^i(*)II^IIII^H 

(3.28) +<Tl(t)\\Vt\\ [ Pt(t ~ s)\\r,(s)\\ d S + Ai Ai (a 1 (t)(-A /l )- 1 T ?t )- 

Jo 

From (|3.18p and from the definition of the non-linear operator, we have 

A 1At ({-A h )- l Vt ) = - b(u hjt , E, (-A h )- l Vt ) - b(u h , E t , (-Ah)- 1 ^) - b(E t , u, (-A^)- 1 ^) 
-ftCE.^C-Ah)- 1 ^) 

= -~{((u h)t • V)E, (-A^)- 1 ^) - ((u M • V)(-A /l )~ 1 ?7 t , E) + (( Uft • V)E t , (-A^ 1 ^) 

-((u fc • V)(-A fe )- 1 r ?i ,E t ) + ((Et • V)u, (-A^ 1 ^) - ((Et • V)(-A^)- 1 r 7t , u) 

(3.29) +((E • V)u tj (-A^)- 1 ^) - ((E • V)(-& ft )~V u*)}. 

For the 2nd, 7th and 8th terms on the right-hand side of (|3.29p , we use Lemmas 13.11 and 13.31 and 
Theorem 3.2 from [8]. For some e > 0, we now find 

^a(t){((u fc>t • V)(-A h )- 1 r ?t , E) - ((E • V)u t , (-A h )- l r, t ) + ((E • V)(-A ft )- 1 ^, u t )} 

(3.30) < C.riCOIlESllllnJlVail^lI^CIIii^tld + 1(11*11!) < «ri(*)ll»I*ll a + C^o-COH^H?.! + iiTe^llEll 2 . 

Similarly, writing E = £ + r) and using Lemmas 13. II and 13.21 we estimate the 4th, 5th and 6th terms 
on the right-hand side of (|3,29p , 

^ a (t) {((u ft • V)(-A /l )- 1 ^,E t ) - ((E* • V)u, (-Aft)- 1 !?*) + ((Et • VX-Aft)- 1 ^,")} 

< CaimttW + WvMmW^ H\\ 1/2 (\H\\i + l|u||i) 

(3.31) < £ffl (t)||77tl| 2 + ^(t)h t ||?-i + ^i(i)ll^l| 2 . 
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To bound the 1st term on the right-hand side of (|3.29p . we first rewrite it as follows: 
(with the notations D{ = and v = (v\,V2)) 

2 r 

((u M - V)E,(-A h )- 1 r ?t ) = V / u hu D i E J ((-K h )- l r lt ) j dx 

2 2 

= -J2 D i u h,t,i E i((-^h)~ l Vt)j dx-J2 / WAAiAKC-Afc)- 1 ^);}^ 

(3.32) = -((V • u M )E, (-A^)- 1 ^) - ((u M • V)(-A ft )- 1 r 7t , E). 
As in (|3.30p . we now have the following bound. 

(3.33) -i CTl (t)((u M • V)E, (-A^)- 1 ^) < evi(t)\\Vt\\ 2 + Ca(t)\\ Vt \\ll + Ke 2at \\E\\ 2 - 
Using similar technique and as in (|3.31j) . we can bound the 3rd term on the right-hand side of (|3.29p . 



1 



(3.34) -- ai (t)((u h ■ V)E t , (-A^)- 1 ^) < eo-i(t)||f| t |I 2 + Ka(t)\\ Vt \\ 2 - 1 + 



|2 



Incorporate (|3.30p - (|3.34p in (|3.29p and then in f|3.28[> . Use Young's inequality and appropriate e > 
and finally integrate with respect time to find 



+ ai(a)R(8)|rcfa<iir / a(s)\\ Vs (s)\\Uds 
Jo Jo 

+K f e 2as {\\<n(s)f + \\E(s)\\ 2 } ds + K f a l (s)U s (s)\?ds. 
Jo Jo 

Use Lemma 5.4 from [8j and (|3.26p to obtain 

(3.35) nmVttl + H f ai(s)\\il s (s)\\ 2 ds<K(t)e 2oct h A + K f a(s)\\ Vs (s)\\-ids. 

Jo Jo 

To estimate the last term of (|3.35|) . we recall the equation in 77 (see (5.12) of [8]): 

(3.36) (rj t , d> h ) + [M(n, <j> h ) + I P(t - s)a(r](s),ct) h ) ds = A ljh (d> h ), 

Jo 

where Ai.^ is given by (I3.16p . We rewrite it using integration by parts. 

(3.37) (r} t ,4> h ) + fia(r),(f> h ) + -fa(fi,(f> h ) + [ &(* - s)a(ri(s),(f) h ) ds = k x , h (4> h ), 

Jo 

where 

fi(t) = j rj(s) ds. 
Jo 

Choose cj> h = cr(t)(- A ft )- 1 T7 i in (13371) . 

^t)\\m\\-i + £ j t W(t)\\v\\ 2 } < f <T t (t)\\n\\ 2 + c*(t)\\rhU{W\i + f p(t - s )||f)( s )||i ds} 

J 

(3.38) +a(t)A hh ((-A h )- 1 rj t ). 
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We estimate the non-linear term Ai^, following the technique used to estimate f|3.29|) . 

a(i)A 1)/l ((-A^)- 1 » 7t ) = -a(t)b(E, u h , (-A^)" 1 ^) - a(t)b(u, E, (-A/j) -1 ^) 

< taj^llEllll^H-idluftlla + ||u|| 2 } 

Incorporate this in (|3.38|) and use kickback argument. We bound and u using Lemma 13.21 and 
Theorem 3.1 from [8], and then integrate with respect to time. 

(3.39) fa^WrjMW-ids + amrjfKc l\ 2as {\\ V (s)\\ 2 + Ms)^ + K\\E(s)\\ 2 } ds. 

Jo Jo 

Apply Lemmas 5.4 and 5.5 from [8] to bound the left-hand side of (|3.39p by i^(t)/i 4 . Incorporate 
this in ()3.35p to complete the rest of the proof. □ 



Lemma 3.5. Under the assumptions of Theorem \3.4\ the following result holds: 

(3.40) e~ 2at I e 2as \\{n-n h ){s)\\lds + e- 2lxt I cri( S )||(u s - M h ,s){s)\\\ds < K(t)h 2 . 

Jo Jo 

Proof. Note that u — = E = £ — 77. And so, keeping in mind (|3.22p and (|3.25p . we only need 
estimates of rj. For that, we choose cj) h = e 2at r] in (|3.36p . 

(3.41) -£{ e 2at \\ v \\ 2 } + fie 2at \\r)\\ 2 + [ /3(t - s)a(r 1 (s),e 2at r 1 ) ds = ae 2at \\r]\\ 2 + e 2at A hh (<n). 
2 dt Jq 

The non-linear term is estimated, similar to (|3.31|) . 

(3.42) A 1A ( V ) = -b(£ + V ,u h , V )-b(u,tv)<4v\\l + K(M\\i + \\v\\ 2 ), 

for some e > 0. Incorporate (|3.17p in (I3.4ip and integrate. Drop the double integral term, following 
Lemma 12.11 

e 2at \\r)\\ 2 + f e 2as \\r){s)\\lds < K t e 2as \\r](s)\\ 2 ds + K f e 2as \\£{s)\\lds. 
Jo Jo Jo 

Use (|3.22p and then apply GronwalPs Lemma to find 

(3.43) \\r}\\ 2 + e~ 2at I e 2Q!S | [ T7 (s) 1 1 ?cZs < K(t)h 2 . 

Jo 

This along with f|3.22|) completes the first part of (|3.40p . 
Next, we choose cj> h = ai(t)r] t in (|3.27p to obtain as in ([3.28H 

lj t {°i(t)\\Vt\\ 2 } + »°i(t)\\Vt\\l <C<rmVt\\ 2 + C<7i{t)\\vh\\Vt\\i 

(3.44) +<n(f)\\Vth / Pt(t ~ s)\\v(s)\\i ds + A 1At (ax(t)ri t ). 

Jo 

For the non-linear term Ai ^j, we recall (|3.29|) . except now, (— Aj l )^ 1 r] t is replaced by rj t . Terms not 
involving Et can be bounded as: (using Theorem 13. ip 



< C(||u M ||i + ||ut||i)||E||i||»7 t ||i < A'(t)(r*(t))- 1 / 2 / l (||u M || 1 + ||u t ||i)||rfc||i. 
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And for terms involving Et, we rewrite it as £ t + r\ t and bound those terms as 

< a|K||i(||€ t ||i + IkllOlkini^lli 72 + ciiuftiiaii^Hidi^lli + htll 1/2 htllf) 

+C(Mth + lll*l| 1/a |l»ltl|} /a )l|u||i||il t || 1/a ||fl t ||; /a + C(||^||i + H^ini^ll^llr/JiHull! 

< ^(ll^llihtlli + HWHWi + IMI 1/2 Nlli /2 )- 

Use Lemma [3.21 and Lemma 3.4 from [8] to bound ||u/j,||i and ||u||i. Now, incorporate the above two 
estimates in Ai^tiVt) an d then in (|3.44[) . Use Young's inequality and kickback argument to find 

^{ CTl (i)||r? t || 2 } + ^xC^H^Hf < C^H^H 2 + C CTl (t)||77||? + iTCt)/. 2 
+Ka 1 (t)(U t \\l + \\ Vt f) + K(t)h 2 a(t)(\M\i + Hill)- 
Integrate with respect to time to get 

<Ti{t)\\mt + ^ f ^WVs^Wids < K f a(s)\\ Vs (s)\\ 2 ds + C f <7i(a)(||f/( a )||? + ds 
Jo Jo Jo 

(3.45) +K(t)h 2 + K(t)h 2 [ <7(s)(||u ft ,,(s)||? + +||u,(s)||?) ds 

Jo 

Use Lemma 13.31 and Theorem 3.2 from [8j to bound the last the term of (I3.45p . Next, we use ()3.25|> 
and (|535P to infer from ([315D 

(3.46) <Ti(i)||»/ t || 2 + /z I a x {s)\\n s {s)\\\ds<K f a(s)\\r, s {s)\\ 2 ds + K{t)h 2 . 

Jo Jo 

To estimate the term on the right-hand side of (|3,46p . we put 4>h = a (t) r lt m (|3-36p . 

(3.47) a(t)\\ Vt \\ 2 + f | WONI?} = ^t(t)NI? " f Pit - s)a( V (s), Vt ) ds + a(t)A lih (rj t ). 
We first rewrite the integral term as: 



o 



a(t) [ Pit- s)a(r](s),'n t ) ds = a(rj t ,a(t) [ /3(t - s)r](s) ds) 
Jo Jo 

= j t [a( V , a(t) j f3(t - s) V (s) ds)} - a(f], <r t (t) J Pit - s) V (s) ds) 
(3.48) -a(ri,<r(t)P(p)Ti) -a(rj,a(t) I /3 t (t - s)rj(s) ds). 



o 



We estimate the non-linear term as (following the idea (|3.32p to avoid 

A llfc (^)<Jr[|f| t [|||E|| 1 (||u fc ||;/ a + ||u||f ). 

Incorporate this along with (|3.48p in (|3.47p . Use Young's inequality and kickback argument. Integrate 
with respect to time to find 

(3.49) a(t)\\v\\l+ fvWWn.WfdsKC f e 2as \\r,(s)\\ 2 ds + K f e 2as \\B(s)\\ 2 ds. 

Jo Jo Jo 

Use the first part of (pTlOj) and (j3"lE3|) to bound the left-hand side of (FTl9"|) by K(t)h 2 . Incorporate 
this in (I3.49p and then in (|3.46p to complete the rest of the proof. □ 
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4 Two-Level Finite Element Method 



In this section, we work with another space discretizing parameter H, that corresponds to a coarse 
mesh. In other words, < h < H and both h, H tend to 0. We introduce associated conforming 
finite element spaces (H#,L#) and (H^,!^) such that (H#,L#) C (H/^L/j). And this two-level 
finite element is to find a pair (u h ,p h ) as follows: 

First Level: We compute the mixed finite element approximation {mhiPh) G (H#,L#) of (u,p) 
of (|2.ip . In other words, we solve the nonlinear problem on a coarse mesh. Find (uff,p#) € (Hjj,ph) 
satisfying 

(u Ht ,4> H ) + fia(u H ,<i>H)+K u H,u H ,4> H ) + [ P(t- s)a(u H (s),4> H ) ds - (p H , V- 0#) = (f,0#), 

•/ o 

(4.1) (V-uj,xh) = 0, 

for (4> h ,xh) e (H. H ,L H ). 

Second Level: We solve a linearized problem on a fine mesh. Find (u /l ,p' 1 ) € (H^, L^) satisfying 

( u ?> 0ft) + 0a(u ft , h )+&(uff, u H , ft ) + / /S(*-a)a(u h (s),0 A ) ds - (p h ,V ■ (j> h ) = (f,0 fc ), 

■/ o 

(4.2) (V-u\ X/l ) = 0, 
for {4>hiXh) e (H ft ,L ft ). 

An equivalent way is look for solution in a weekly divergent free space. 
First Level: Find u# 6 J# satisfying 

(4.3) (u Ht ,cf) H ) + na(u H ,(j) H ) + b(u H ,u H ,cf) H ) + / fi(t - s)a(u H (s),<j> H ) ds = (i,(f> H ), 

Jo 

for 4> H £3 H . 

Second Level: With ujj being the solution of (14. 3D . find u' 1 E satisfying 

(4.4) (u t \ h ) + /ia(u ft , 4> h ) + 6(u Hl u H , fc ) + f P(t - s)a{u h (s),cb h ) ds = (f, h ), 

■/ o 

for <p h €J h . 



Remark 4.1. XTie well-posedness of the above systems can be seen from the facts that ( ETTp or 
is Galerkin finite element approximation and hence is well-posed as is stated in Section 3. And ft4-2\ ) 
or ^4-4\ ) represent linearized version. Given uh and with suitable u h (0), it is therefore an easy task 
to show the existence of an unique solution pair (u h ,p h ) (or an unique solution u h ) for the linearized 
problem following the foot-steps of the non-linear problem. 

Following Lemma 13.21 we can easily obtain the a priori estimates of u#. 

Lemma 4.1. Under the assumptions of Lemma \3.S\ and for u#(0) = Phuq, the solution vlh of |j5p 
satisfies, for t > 0, 

(4.5) \\u H (t)\\ + e~ 2at f e 2as \\u H (t)\\ 2 ds < K, 

Jo 

(4.6) \\u H (t)\\ 1 + e- 2at f e 2as \\u H (t)g ds < K, 

Jo 

(4-7) (T*(t)y/ 2 \\u H (t)\\ 2 <K, 

where r*(t) = min{l,t} and K depends only on the given data. In particular, K is independent of 
H and t. 
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The following higher-order estimate of u# is required for error analysis. The proof of the same is 
similar to that of Lemma 13.31 



Lemma 4.2. Under the assumptions of Lemma \4.1\ the solution u# of |J~5p satisfies, fort > 0, 
(4.8) (r*(t)) 2 \\u Hit (t)\\ 2 + (r*(t)) 3 ||u^ t (i)||2 + e ~ 2at ? » 2as <-*> 



f e 2as (T*{s))\\^ s (s)\\\ds<K. 
Jo 



Below, we deal with the a priori estimates of u . 

Lemma 4.3. Under the assumptions of Lemma \4.1\ the solution u h of lji4-4\ ) satisfies, for t > 0, 



(4.9) 



\u h (t)\\+e- 2at I' e 2as \\u h {t)\\ 2 ds < K, 
Jo 



(4.10) ||u ft (t)||i + e" 2 "* f e 2as \\u h (t)\\ 2 ds < K, 

Jo 

Proof. Given u# with Lemma (|4.ip . we choose 4> h = e 2at u h (t) = e at u h (t) in ([4.4p to obtain 
^||iiY-a[|u h || 2 + /x||u h ||?+7 /"e-^^aCu^W.u") ds 

— ■ (XL 







(4.11) =(f,u h )-e 2at b(u H ,u H ,u h ) 

Use Cauchy-Schwarz inequality, Poincare inequality with first eigenvalue of Stokes operator as the 
constant and Young's inequality to have 

(4-12) (f,u h ) < ||f||||u /l || < -^llflHIu^lli < ^\\u h \\l + ^||f|| 2 . 

X{ f 4 

From Lemma 13.11 and Lemma I4.1|. we find that 

,/ h\ ^ II Ml/2 II 1 1 3/2 r . M l/2|| 1 1 1/2 

b(ufl-,UH,u*) < ||uff|| 1 ||uh||i II u/x. II 1 \\VL h \\{ + IIu^IIIIuhIIiIIu^Hi 



<K + ^\\u h \\ 2 
4" 



Therefore, 



(4.13) e 2a *&(uff,u/r,u h ) < Ke 2at + j\\n h \\\. 
Putting the estimates (|4.12|) - ()4.13|) in (|4.1ip gives us 

-||u h || 2 - 2a\\u h \\ 2 + fiWuX + 2 7 [* e^ 5 - a ^a(u h (s),u h ) ds < Ke 2at + — ||f|| 2 . 
at Jo /uAi 

Use Poincare inequality to get 

(4.14) ^\\u h \\ 2 + U- ^)\\n h \\{ + 2 7 f 'e-( s - a ^a(u h (s),u h ) ds < Ke 2at + -^||f|| 2 . 
at Ai Jo /xAi 

Since < a < min{^Ai/2, <5}, we have that fi — 2a/ Xi > 0. 

Integrate (|4.14p with respect to time and use the positivity of the kernel f3, that is, Lemma 12.11 to 
drop the resulting double integral. This results in 

HuW + t/x-^) /Vwilf^lCe 2 "*. 
A i Jo 
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Multiply by e~ 2at to conclude (jlH . 

For the next estimate, we choose cj) h = e 2at AhU h (t) in (|4.4p and proceed as above. For the non-linear 
term, by the definition of the operator b, we have 

b(u H ,u H , A h u h ) = -((u H ■ V)u H ,A h u h ) - -((u H • V)A h u h ,u H ). 
Use Lemma 13. II to obtain 

(4-15) ^((u H -V)u H ,A h u h ) <C||u h || 1/2 ||uh||2 /2 ||u^||i||u /i || 2 . 

Following the idea (|3.32p . we have 

((Ufl • V)A A u h , u H ) = -((V • u H )u H , A h u h ) - (( Uff • V)u H , A h u h ) 
(4.16) <C||uh|| 1/2 ||uh||2 /2 ||uh||i||u' 1 || 2 . 

These estimates let us have, with cj) h = e 2at A/ l u /l (t) in (|4.4p : 

u h \\l - a\\u h \\ 2 + fi\\u h \\ 2 + 7 / e^ 5 - a ^- s \A h u h (s), A h u h ) ds 

Jo 

(4-17) < ||f||||u h ||2 + Ce a *||u ff || 1 / 2 ||u i/ ||^ /2 ||u ff ||i||u' l ||2. 

Using Young's inequality and Poincare inequality as earlier, we find that 

+ ^)||^||l-h2^^*e-^-«)(*--)(A ft ^( S ), A^ ft ) ds 

(4-18) < C||f || 2 + ce 2at {\\u H \\ 2 \\u H \\i + \\u H \\ 2 }. 

Integrate (|4.18p with respect to time, use Lemma l4.1l and finally multiply by e~ 2at to obtain (|4,10p , 

□ 

5 Error Estimate 

In this section, we present the error estimate for the spatial approximation, that is, two-level finite 
element Galerkin approximation. We achieve our desired results through a series of Lemmas. We 
denote the error, due to two-level method, as e = — u . 
From the equations (|3.3p and (|4.4p . we have the following error equation: 

(5.1) (et, h ) + /ia(e, fc ) + / P(t - s)a(e{s),(f> h ) ds = A h {<f> h ) V0 ft G J h , 

Jo 

where 

^h((f>h) = b(u H ,u H , 4> h ) - b(u h , u h , 4> h ) 

(5.2) = b(u H , u H - u h , 4> h ) + b(u H - u h , u h , <j> h ). 

Lemma 5.1. Under the assumptions of Theorem \3. 1\ and with the additional assumption that u ft (0) = 
Ufe(O) and for < a < min{//Ai/2, 5}, the following estimate 

(5.3) l|e(t)|| 2 i + e~ 2at f e 2as \\e{s)\\ 2 ds < K(t)H 4 

Jo 

holds, for t > 0. 
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Proof. Choose <f> h = e 2at (-Aft) _1 e(t) = e at {- h h )~ l e(t) in to obtain 
(5.4) I^||e|| 2 _ 1 -a||e|| 2 _ 1 + / ,||e|| 2 + 7 f e'^^ (e( S ), e) = e^tf-A,)-^) 



2 j i0 
From (|5.2p and Lemmas 13.11 13.21 and 14.11 we have 

A((-A /l )" 1 e) = b{u H , u H - Ufc, (-A A ) -1 e) + 6(u^ - u h , U/l , (-A^e) 
< C\\u H - u ft ||||e||{||ui/||i + llu/Jli} < K\\u H - u h \\\\e\\. 

We have again used the idea (|3.32j> . Incorporate this in (|5.4p . Use Young's inequality and kickback 
argument. 

rt 



(5.5) ^-||e|| 2 _ 1 + /U ||e|| 2 + 2 7 / e"^^) (e(s), e) ds < 2a||e|| 2 _ 1 + K\\u H - u h f 

at J k) 



Integrate (I5.5P with respect to time. The resulting double integral is dropped following Lemma l2.ll 
Observe that 

(5.6) f \\(u H -u h )(s)\\ 2 ds< [ t \\(u-u H )( s )fds+ f \\(u-u h )(s)\\ 2 ds<Ke 2at H 4 . 
Jo Jo Jo 

We have used here the Lemma 5.4 of [8]. And hence 

(5.7) ||e(t)||_i + n I ||e(s)|| 2 ds < 2a f \\e(s)\\-ids + Ke 2at H A . 

Jo Jo 

Multiply by e~ 2at and use Gronwall's Lemma to complete the rest of the proof. □ 

Remark 5.1. Due to non-smooth initial data, in the estimate of ||e(t)||, it is essential to introduce 
a(t) = T*(t)e 2at , r*(t) = min{l,t} ; in order to avoid nonlocal compatibility conditions (see J2|/ and 
references cited there). However, the resulting integral term (with cb h = a(t)e(t)) is now no longer 
non-negative after integration. A direct estimate of this term is not possible. To handle this new 
difficulty, we present the following notation 



e(t) = / e(s) ds 
Jo 



and obtain an estimate for 

rt 

O 2o:i||~/„M|2, 



e 2at \\e{s)\\(ds. 

And with this intermediate estimate, we manage to estimate the integral term with no loss of rate of 
convergence. 

First, we integrate the error equation (|5.ip in time and observe that Gt(t) = e(t), since e(0) = 0. 
Also from (5.18) of [SJ, we find that the resulting double integral can be written in single integral 
form. 

(5.8) (&t,4> h ) + fia(e,4> h )+ [ P(t - s)a(e(s), <f> h ) ds = f A h (d> h ) ds V0 fo € 3 h , 

Jo Jo 

Lemma 5.2. Under the assumptions of Lemma \5.1\ the following estimate 

(5.9) ||e(t)|| 2 + e- 2a ' f e 2as \\e{s)\\\ds < Ke 2at t l/2 H A 

Jo 

holds, for t > 0. 
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Proof. Put (f> h = e 2at e(t) in §B) to find 

I^.| e 2«t||g||2| _ ae 2at||e||2 + Me 2^||g||2 + ^ f* e -0-a)(*-«) o ( e <"g( s ) ) e <rfg) ds 

2 at Jo 
(5.10) = e 2crf / A ft (e) ds. 



Now from (|5.2p and Lemma 13.11 we have (following (|3.32p ) 

Afe(e) = 6(u#, u H - u h ,e) + 6(u# - u h , u h ,e) 
( 5 - n ) < C\\u H - Uftllllellilllu^llJ^lluj/H^ 2 + llufellJ^Hufell^ 2 }. 

Incorporate (|5.1ip to obtain from (|5.10p 

-^{e 2a *||e|| 2 } -2ae 2Q *||e|| 2 + 2/ie 2a *||e|| 2 + 27 f e^ 5 -^^ a(e as e(s), e at e) ds 



dt 



< ee^HelU / ||(u* - u h )(s)\\{\\u H ( S )\\l /2 \\u H (s)\\l /2 + \\u h (s)\\l /2 ||u,( S )||f } ds. 
J o 

Use Young's inequality, kickback argument and then Poincare inequality. Again apply Young's 
inequality and then Lemmas 14.11 and 13.21 to find 

^{ e 2at||ef} + ( // _^) e 2at||g||2 + 27 r e ~(<5-«)(t-) a(e -e( S ),e Q 'e) ds 
dt ^i J t Q 

<Ce 2at \ r||(u fl -u fc )( a )||{||u ff ( S )||l /2 ||u ff ( S )||5 /2 + ||u fc ( S )||} /2 \\Ms)\\T } ds}" 
Jo J 

<Ke 2at f \\{vL H -u h )(s)\\ 2 ds f \\{\\uH{s)h + \\u h {s)\\ 2 } ds 
Jo Jo 

(5.12) <Ke 2at flKuH-u^fds f s~ 1/2 ds < Ke 2at t 1/2 f \\(u H - u h )(s)\\ 2 ds. 

Jo Jo Jo 

We use the fact that 1 < e at since a and t are positive. Use (|5.6p to get 

f ||(utf - u h )(s)\\ 2 ds < f \\(u H - u h )(s)\\ 2 ds < Ke 2at H\ 
Jo Jo 

With this, we now integrate (|5.12p in time and drop the double integral, since it is non-negative, see 
Lemma 12.11 

(5.13) e 2«*||g( t )||a + (^_^) f e 2as \\e(s)\\ 2 ds<Ke 4at t l l 2 H\ 

^i Jo 

Multiply (|5.13p by e~ 2at to complete the rest of the proof. □ 
Lemma 5.3. Under the assumptions of Lemma \5.1\ the following estimate 



(5.14) r*(t)||e(t)|| 2 + e~ 2at f a(s)\\e(s)\\ 2 ds < K{t)H* 

Jo 

holds, for t > 0, where a(t) = T*(t)e 2at , r*(i) = min{l,t}. 
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Proof. We first recall that f3(t) = je st and now use integration by parts to find 
(5.15) / P(t - s)e(s) ds = je(t) -5 [ f3{t - s)e(s) ds. 



Keeping this in mind, we choose 4> h = a(t)e(t) in (|5.ip to obtain 

^{a(t)||e|| 2 } + 2/i a(t)\\e\\l = a t (t)\\e\\ 2 + 2<r(t)A ft (e) - 2 7 a(t)a(e,e) 

(5.16) +28a(t) f f3(t- s)a(e(s),e) ds. 

Jo 

We repeat the non-linear estimate (15. lip . 

Afc(e) < C\\u H - u^llllellilllu^llJ^HuHll^ 2 + ||ufc||i /2 ||ufc||2 }• 
And hence, we find from (|5.16p 

-{<7(t)||e|| 2 } + 2fi a(t)\\e\\ 2 < (1 + 2a)||e|| 2 + 2<r(t)||e||i{7||e||i + 5 J (3(t - s)||e||ids} 

+Ccr(i)||e||i||uH - UftlKlluij-ll^llujj-H^ 2 + ||uh||J' 2 ||u h ||2 /2 }. 
Use Young's inequality and then kickback argument to yield 

^{a(i)||e|| 2 } + ft a(t)\\e\\ 2 < (1 + 2a)||e|| 2 + Ca(t){ 7 ||S||i + S f*P(t ~ s )Mi ds Y 

(5.17) +Ca{t)\\u H - u h \\ 2 {\\u H \\i\\u H \\ 2 + \\u h \\i\\u h \\ 2 } . 
Use Lemmas 14.11 and 13.21 to estimate the last term of (|5.17p as follows: 

Ca(t)\\u H - u/j|j 2 {||u//||i||u//|| 2 + ||u A ||i||u A || 2 } < K\\u H - u h \\ 2 . 

We have used the fact that r*(i) < 1. Incorporate this in (|5.17p and integrate in time. Repeat the 
technique of (|5.6|) . 

2 ... I _/„MI„/'„M|2 t / n I J ||A/„\||2 , „2as\\~f„\\\2\ j„ , V „2at ttA 



a{s)\\e(s)\\\ds < C J* {||e(s)|| 2 + e 2QS ||e(s)|| 2 } ds + Ke 2at H 4 



(5.18) o-(t)\\e(t)\\ 2 + n / a(s)\\e(s)\\(ds<C 

Jo 

Now use Lemmas 15.11 and 15.21 and then multiply the resulting inequality by e~ 2at to complete the 
rest of the proof. □ 

Lemma 5.4. Under the assumptions of Lemma \5.3\ the following holds: 

(5.19) (r^t^Wemi^mH 2 . 
Proof. Put 4> h = o-(t)e t (t) = er(i)e(t) in ((OJ to find 

(5.20) 

^)l|e*|| 2 + f|{^WII^II 2 } " f "tWIN? = fft* ~ «)«(S(-),e) ds + a{t) j\ h (e) ds. 

As in (|5.1ip . we have, using Lemmas 14.11 and 13.21 

(5.21) a{t) f Afc(e) ds < Ka(t)\\e\\i I \\u H {s) - u h (s)\\{ ||u H (s)||^ /2 + K(s)|| 2 /2 } ds. 
Jo Jo 
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Put (|5.2ip in (|5.20p and integrate with respect to time. Use Cauchy-Schwarz and Young's inequalities. 
For the resulting double integral term proceed as in (|5.12p and f|5. 13|) to observe that 

(5.22) M*)l|e||i + 2 f a(s)\\e s (s)\\ 2 ds <c f e 2as \\e{s)\\\ds + c / a{s) ||e||fds + Ke 4at t 3 / 2 H 4 . 

Jo Jo Jo 

Use Lemmas 15.21 and 15.31 to complete the rest of the proof. □ 
Lemma 5.5. Under the assumptions of Lemma \5.1\ the following estimate 

(5.23) (r*(t)) 2 \\e(t)\\l + e- 2at ai(s)\\e s (s)\\ 2 ds < K{t)H A 

Jo 

holds, fort>0, where a^t) = (t* (t)) 2 e 2at . 

Proof. Choose 4> h = oi{t)e t {t) in (|5.1|) and incorporate (|5.15|) to obtain 

^iWlletll 2 + | ^ W*)||e||?} = f <xi, t (i)||e||? + ffi(t)A h (et) - 7 <ri(i)a(e,et) 

(5.24) +50-^)1 (3(t- s)a(e(s),e t ) ds. 



o 



Multiply (|5.24p by 2 and integrate with respect to time. Use integration by parts to observe that 

- 2 7 / a 1 (s)a(e(s),e s ) ds 
Jo 

= -2 7 (Ti(t)a(e,e) + 2 7 / a ljS {s)a(e(s), e) ds + 2 7 / <ri(a)||e(a)||? ds 

Jo Jo 



(5.25) <T^i(t)IH|? + a7i(t)||e||f + C / a(s)\\e(s)\\(ds + C a(s)\\e(s)\\(ds. 



4 



Jo 



And 



25 / (Ti(s) / (3(s — r)a(e(r),e s ) c?t <is 
io io 



=2<5cri (i) / /?(t - s)a(e(s),e) ds - 25 / cr liS (s) / - r)a(e(r), e(s) dr ds 
Jo Jo Jo 

+ 2<5 2 / ^(s) / /3(s -r)a(e(r),e(s)) dr ds - 2<5 7 / <r 1 (s)a(e(s), e(s)) 
io io io 

(5.26) ^^(tJHeHf + c/ e 2as \\e(s)\\ 2 ds + C [ a(s)\\e(s)\\ 2 ds. 

4 io io 

We have again used the fact that r*(t) < 1. And now, after integration of (|5.24|) . we incorporate 
([Q5> (15361) to find 

2 / (7i(*)K( fl )|| 2 tfa + ^(7i(t)||e||?<C f e 2as \\~e{s)\\ 2 ds + C f a(s)\\e(s)\\ 2 ds 
Jo 3 Jo io 

(5.27) + C<7 1 (i)||e||? + 2 / u x {s)K h {e s ) ds. 

Jo 

We rewrite the non-linear terms as follows: 

A h (e t ) = b(u H , u H - u h ,e t ) + b(u H - u h , u h , e t ) 

= ^{b(u H ,u H - u h ,e) + b(u H - u^u^e)} - b(u Htt ,u H - u h ,e) 
- b(u H , u ff i - Uh,t, e) - b(u H ,t - ^h,u u h, e) - b(u H - u h , u htt , e). 
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And hence 

cr 1 (t)A ft (e t ) = ai(t)(b(u H ,u H - u h ,e t ) + b(u H - u h ,u h ,e t )) 

= ^{cri(t)(6(u H) UH - u h ,e) + 6(uh- - u h ,u h ,e))} 

- ax,t(t)(b(uH,UH - u/t,e) + 6(utf - u/^u^e)) - cri(t)6(uH ( t,UH - u h ,e) 

(5.28) - ffi(t)(6(u ff ,Uff i( - u^t,e) + b(u H ,t - u h;t ,u h ,e) + b(u H - u h ,u h;t ,e)). 

As seen in (|5.11j) . we have 

- b(u H) t, uh - u ft , e) - b(u H - u h , u h>t , e) 
<C||uh - Ufellllellilllu^illi^Hu^tlla' 2 + ||u^,t||i Hu^tH^ 2 } 

Use Lemmas 13.31 and 14.21 to conclude that 

(5.29) - b(u Hit ,u H - u h ,e) - b(u H - u h ,u hjt ,e) < K(T*(t))~ 3/2 \\u H - u^||||e||i. 
Similarly 

(5.30) - b(u H ,u H - u h ,e) - b(u H - u h ,u h ,e) < K(t* (t)y 1/2 \\u H - u h \\\\e\\i 
and 

(5.31) - b(u H ,u H ,t - u hit ,e) - b(u H ,t - u htt ,u h ,e) < K(r*(t))~ 1/2 \\u H ,t - u h ,t\\ ||e||i. 
Incorporate (|5.29p - (|5.3ip in (|5.28j) and integrate with respect to time. Re-use (15.30P to find 

/ ai(s)A h (e s ) ds <a t (t) (b(u H , u H - u h ,e) + b(u H - u h ,u h ,e)) 
Jo 

+K f e 2ccs {(T*{s)) 1/2 \\u H - u h \\ + (r*( S )) 3 / 2 ||u^ -u M ||}||e( S )|| 1( i S 
J o 

<Ke 2at (T*(t)f/ 2 \\u H - u h \\\\e\\i + C f a(s)\\e(s)\\ 2 ds 

Jo 

+ K [\ 2as {\\u H -u h \\ 2 + (T*(s)) 2 \\u H , t -u h , t \\ 2 } ds 
Jo 

<^<Ti(t)\Wi+K<r(t)\\vH-vi h \\ 2 + C / a( S )||e( S )|| 2 d S 
6 Jo 

(5.32) +K [ e 2as {\\u H - u h \\ 2 + (r*( S )) 2 ||u^ t - u M || 2 } ds. 

Jo 

Put ([5321 in (pHTD to obtain 

(5.33) 2 raiW||e,( 8 )|| 2 d S + ^<7 1 (*)||e||? <C f e 2as \\e(s)\\ 2 ds + C I a(s)\\e(s)\\ 2 ds 

Jo 3 Jo Jo 

+C<7 1 (t)||e|| 2 + Ka(t)\\u H - u h \\ 2 + K [ e 2as {\\u H - u h \\ 2 + (r*( S )) 2 ||u Hjt - u M || 2 } ds. 

Jo 

Now, use Lemmas 15.21 15.31 15.41 and 13.41 For the remaining part, we use (|5.6p and Lemma 5.4 and 
Theorem 5.1 from [8]. This completes the rest of the proof. □ 

We present below, the pressure error estimate. 
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Lemma 5.6. Under the assumptions of Lemma \5.1\ and additionally that the assumption (B2') 
holds, we have 

(5.34) llfa-ftmKKWT*®)- 1 !?. 
Proof. The LBB condition (B2') tells us that, for t > 

, K «« II fen ^ n (p h -p h ,V ■ (f> h ) 
(5-35) \\p h -p \\ L 2/ R < C sup jj— r . 

o^ h eH„ \\<Phh 

From (J32J) and (JOJ), we have, for h G H h 

(5.36) (p h - p h , V • h ) = (e t , 4> h ) + na(e, (j> h ) + f /3(t - s)a(e(s), h ) ds - A h (<p h ). 

Jo 

Using (I5.15P and (15. 11 H . we obtain from (I5.36P 

(p h -p\ V • cj> h ) < CH^IIidletll-Lh + Helix + ||e||i + {e~ 2at f e 2as \\e(s)\\lds) 1/2 } 

Jo 

(5.37) +K(T*(t))- 1 /*\\u H - UfcHII^Hi, 
where 

n I, <et,4>h> 
\\et\\-i,h= sup — t— J, . 

o#0 fc eH h II ^ Hi 
Taking supremum over a bigger set, we find that 

(r. qq\ II II / II II < e *' ^ > 

( 5 - 38 ) ll e t||-i,ft < ll e tll-i = SU P — n-rn — • 

o^eHj II0IU 

Incorporate this in (|5.37|) and now, for / (f> h , divide by ||<£/J|i. Use (|5.35|) and then Lemmas 15 .51 
E2Dand[5J2]to get 

(5.39) \\ Ph - P h \\ L 2 /m < C\\e t \U + K{t){r*{t))- l H 2 . 



As usual, we have rewritten uh — as (u — u/J — (u — u^) and used Theorem 13.11 
For the negative norm, we have, for <fi € Hj 

(e t ,0) = (et,<f>-P h <f>) + (et,P h <f>). 

Since £ J/,, we use the error equation f)5. 1 j) 



(e t , 0) = (e t , - - pa(e, P h 0) - f p(t - s)a(e(s), P h (f>) ds + A h (P h </)), 

Jo 

where is given by (|5.2p . Using discrete incompressibility condition, ^T 1 -stability of Ph, (Bl) and 
(15. lip , we obtain 

(e t ,0) < C/ l ||e t ||||0|| 1 + C(||e|| 1 + f (3(t - s) \\e(s) hds) ||^||i 

j 

(5-40) +C||uH-u /l ||(||u // ||l /2 ||u^||^ /2 + ||u // ||l /2 ||u ff ||^ /2 )||0|| 1 
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We can rewrite the integral term as in (|5.37p and use the estimates used in (|5.37p to obtain from 

||et||_i <Ch\\e t \\+K(t)(r*{t)y l H 2 . 

Incorporate this in (|5.39p . 

(5.41) \\ Ph - p h \\ L 2/^ < CH\\e t \\ + K(t)(r*(t))- l H 2 . 
Assuming the following sub-optimal estimate (since it is of O(H)) 

(5.42) ||et|| < ^(r^t))- 1 ^, 

we complete the rest of the proof. □ 

Remark 5.2. It is sufficient to have a sub-optimal estimate of \\et\\ to attain our result. Instead, 
had we chosen to estimate the first term on the right-hand side of A5.36\) as 

(et,<t>h) < l|et||||0ftll < C||e t ||||^||i, 

we would require an optimal estimate of \\et\\, which would be 0((r*(i))~ 3 / 2 if 2 ). And then the 
pressure error estimate would read 0((t* (t))~ 3 ^ 2 H 2 ) , which is worse that the estimate we have 
attained here, that is [5.3$ - 

Lemma 5.7. Under the assumptions of Lemma \5.1\ the following holds: 

(5.43) ||et|| ^K^it))- 1 !!. 

Proof. We differentiate the error equation (15. ip with respect to time and then incorporate (|5.15p . 

(ett, <f> h ) + na(et, <f> h ) + -ya(e, d> h ) - Sja(e, <f> h ) + 5 2 f3(t - s)a(e(s), (f> h ) ds 

Jo 

(5-44) =A M (0J V^6 4 

where 

J^h,t{4>h) =K u H,t, uh - ™h, 4> h ) + b(u H , u Hit - u hit , 4> h ) + b{u Htt - u htt , u h , 4> h ) 
(5.45) +b{u H -u h ,u htt ,(f> h ). 

Choose 4> h = a\{t)et in (|5.44jl to find 

^{^(OiledH+^iWIIedli^^lletf + cTiWA.^e,) 

(5-46) c<ri(t)||e t ||i{||e||i + ||e||i+ f 0(t - a)||e(*)||ida}. 

Jo 

Using Lemma |3.1| we observe that 

b(u H ,t,UH ~ u-h,cri{t)et) = -<7i(i){((uff ]t • V)(u# - u h ),e t ) - ({u H ,t • V)e t ,u H - u h )} 

< Co-i(t)||uH,t|| 1/2 |[uH,t||i /2 (||uH - u^||i||et|| 1/2 ||e t ||J /2 + ||e t ||i||uij - u h || 1/2 ||u^ - bu h \\{ /2 ) 
(5.47) < Ccri(t)||u.H,i||l||uff - lift. 1 1 1 1 1 e t 1 1 1 . 
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Repeating (|5.47p . we obtain 

&i(t)A h;t < Cai(t)(\\u H ,t\\i + llu^tllx)!!^ - u/JiHetHi 
(5-48) Co"i(*)(||uh||i + ||ufe||i)||uH,i - u^illilletHi 

Incorporate (|5.48p in (|5.46p . Use Young's inequality and kickback argument. Then integrate the 
resulting inequality to observe that 

a^WetW 2 + fi f t *i(s)\\e a (a)\\lda<C f a(s)\\e s (s)\\ 2 ds + C I ffi(s){||e(a)||? + ||e(*)||?} ds 
Jo Jo Jo 

+C f e 2as \\~e(s)\\ 2 ds + K{t)H 2 I a{s){\\u H , s {s)\\l + ||u M ( S )|| 2 ) ds 
Jo Jo 

(5.49) +K (7i(s)||(uH, s -Ufc, 8 )(a)||?£fa. 

Jo 

We have used Lemmas 13.21 and FPU to bound ||u#||i and 1 1 u/i 1 1 1 . As earlier, we have rewritten u^ — 
as (u — Uh) — (u — u#) and have used Theorem 13. 11 Use Lemmas 13.31 and 14.21 to bound ||iiH,t||i and 
||u/i t||i under integral on the right-hand side of (|5,49p . Next, use Lemmas 15.31 and 15.21 Use Lemma 
23 to estimate the last term on the right-hand side of (|5.49ft . 



(5.50) ai(t)||e t || 2 + f a 1 {s)\\B s {s)\\lds < C f ' a{s)\\e s (s)\\ 2 ds + K(t)H 2 . 

Jo Jo 

To estimate the term on the right-hand side of (|5,50p . we put <j> h = a(t)et in (|5.ip . 

(5.51) a(t)\\et\\ 2 + \ ^{<x(i)||e|| 2 } = \a t {t)\\e\\\ - a{t) j* (3(t - s)a(e(s),e t ) ds + a(t)A h (e t ) 
We handle the integral term as in (|3.48p . We take care of the non-linear term as earlier. 

A h (e t ) < K||et||||uH - UftlliXHu/Jla' 2 + llu^ll^ 2 )- 
Incorporate this in (|5.5ip and after simplifying, integrate with respect to time. 

(5.52) v{t)\\e\\l+ f a(s)\\e s {s)\\ 2 ds <C f e 2as \\e{s)\\\ds + K f e 2as \\E{s)\\\ds. 

Jo Jo Jo 

Apply Lemma 13.51 to obtain 

(5.53) *(*)l|e||i + f <r(s)\\e s (s)\\ 2 ds < C f e 2as \\e{s)\\lds + K{t)H 2 . 

Jo Jo 

For the term on the right-hand side of (|5.53p . we put cf) h = e 2at e in (|5.ip to find 



1 d 

^ 2dt 



{e 2a '||ej| 2 } + ne 2at \\e\\i + f (3(t - s)a(e{s),e 2at e) ds = ae 2at \\e\\ 2 + e 2at A h {e). 

Jo 



Note that 

Ah(e) < (\\uhWx + ||uH-||i)||e||i||u H - u/Jli < i^||e||i||u^ - u h \\x. 
We incorporate this in ()5.54p and then integrate with respect to time. Drop the double integral term. 



e 2at 



|e|| 2 + f e 2as \\e{s)\\lds<C f e 2as \\e(s)\\ 2 ds + K [ e 2as \\(u H - U?i )(s)|| 2 ds. 
Jo Jo Jo 



Use Lemmas to observe that 

I ■ r ■ 3 



+ e- 2at r e 2as ||e(s)|| 2 d S < K(t)H 2 
Jo 



Use this to complete the estimate in (I5.53|) . which, we use, in turn, in (15.501) to complete the rest of 
the proof. □ 



23 



Remark 5.3. In the Lemmas fff.^D and \5. 6\) , we have obtained 

(5.55) ||K " u h )(t)||i + \\(Ph ~ P h )(t)\\ < K{t)(r*{t))- l H\ 
Combining this with Theorem \3.1\ we have 

(5.56) ||(u - u ft )(i)||i + \\(p- P h )(t)\\ < K(t){{r*{t))- l l 2 h + (r* (t))' 1 H 2 } . 

We note that for H = 0{h 1 / 2 ), both Galerkin FE and two-level methods produce similar results, away 
from initial time. But while the first method solves the full non-linear equation on a grid of mesh-size 
h, the second only does that on a coarse mesh (grid of mesh- size H, H « h) and supplements it by 
solving a linearized problem on the fine grid (mesh-size h). 

Remark 5.4. Note that the improvement is shown only in energy norm. In other words, in the L 2 - 
norm of velocity error, we see no improvement of the two-level over Galerkin FE. The result achieved 
here, i.e, \5. 55}) . is nothing but super-convergence result and hence, we can not expect improved 
estimate for I? velocity error (at least not by the proof techniques used here). We would like to point 
that out that this super- convergence is also established in jllOf for Navier-Stokes equations, but for 
non-linear Galerkin method, which is similar in nature to two-level or two-grid methods. And in 
f^, a static two-grid method (a discrete steady Oseen-type problem is considered at the second level 
) is discussed for Navier-Stokes equations, where numerical results show that, for mini element, I? 
velocity error has same rate of convergence for both Galerkin FE and the static two-grid methods. 
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